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1 Introduction 



Effective cfiiral quark models liave been used extensively to describe the low- 
energy phenomena associated with the dynamical breaking of the chiral symme- 
try. Of particular interest are models which include the Dirac sea of quarks, since 
they allow for a common description of mesons and baryons. The models share 
the feature of having an attractive four-quark interaction which is either: 

- postulated ||I|]; 

- derived from a model of the underlying QCD structure of the vacuum |0,||,§J^; 

- obtained by modeling the effective gluon propagator and using Schwinger- 
Dyson techniques 0,0; 

- derived from QCD inspired models P|,P| JTI1| , P!T| . 

For extensive reviews with a focus on solitons, and numerous additional ref- 
erences, see 0,0JT^. The four-quark interaction leads to the chiral symmetry 
breaking, which is, in the framework of these models, the key dynamical ingredi- 
ent. 



The various models fall into two categories: local models, where the four-quark 
interaction is point-like and where ultra-violet divergences are removed by in- 
troducing cut-offs in the quark loop, and nonlocal models, where the four-quark 
interaction is smeared, such that all Feynman diagrams in the theory are fi- 
nite. Nonlocal models arise naturally in several approaches to low-energy quark 
dynamics, such as the instanton-liquid model fl^ or Schwinger-Dyson resumma- 
tions [^. For the derivations and farther applications of nonlocal quark models 
see, e.(/.,p,|15i|T6i|7i|8i|19Hg,0,|3,|l,|lj2|j2|j23j2|j2i. 



Considerable effort has been exerted to describe baryons as solitons of effective 

The solitons appear as bound states 



chiral quark models ||T2| , |30| , |31| , |32| , |33| , |34| , |35| , |36| 
of quarks. Except for our earlier- reported work p7| , |38[] , solitons have so far only 
been obtained from local models, which typically use the Schwinger proper-time 
or the Pauli-Villars regularization of the Dirac sea |T^ , |5U[ . One problem encoun- 
tered with the proper-time regularization, by far the most commonly used, is 
that the solitons turn out to be unstable unless the auxiliary sigma and pion 
fields, introduced in the process of semi-bosonization, are constrained to lie on 
the chiral circle [EO|,ETi. Such a constraint is external to the model itself. Nonlo- 



cal models do not suffer from this instability: we have shown p7| , |38| that stable 
solitons exist in a chiral quark model with nonlocal interactions without the ex- 
tra chiral-circle constraint. In fact, we find appreciable deviations from the chiral 
circle. Furthermore, in local models, the regularization is applied only to the real 
(non-anomalous) part of the Euclidean quark-loop term. The finite imaginary 



(anomalous) part is left unregularized in order to properly describe anomalous 
processes [|^. One may view such a division as rather artificial, and we find 
it quite appealing that, with nonlocal regulators, both the real and imaginary 
parts of the action are treated on equal footing. Moreover, the anomalies are 
preserved [p^p!7| , p!9| , p^ , p5[1 and the charges are properly quantized [^. Finally, the 
momentum-dependent regulator makes the theory finite to all orders in the 1/Nc 
expansion. This is in contrast to local models, where inclusion of higher-order- 
loop effects requires extra regulators |4^J4^J45|J46| , ^ and the predictive power is 
diminished. 

The above-mentioned virtues of nonlocal models do not come for free. Two com- 
plications arise. First, soliton calculations require an extra integration over an 
energy variable, which has to be performed numerically. Second, the Noether cur- 
rents in nonlocal models acquire extra contributions, which are, however, needed 
to conserve the Noether currents and anomalies. Nevertheless, an ambiguity re- 
mains: the transverse parts of currents are not fixed by current conservation 
and their choice is effectively part of the model building. This problem has been 
known for a long time in nuclear physics, where the transverse part of the meson- 



exchange currents is ambiguous p0| , |211| . It is not possible to get rid of this am- 
biguity in gauging nonlocal models. An ideal solution would be to first gauge 
the underlying theory {e.g. the instanton model of the QCD vacuum), and then 
to derive an effective gauged quark model. This has not been attempted so far, 
so that in effective nonlocal quark models we need to deal with transverse cur- 
rents which are not uniquely defined. Fixing these currents is a part of the model 
building. 

In this paper we demonstrate the existence of solitons and we give the description 
of formal and technical aspects of constructing solitons in nonlocal chiral quark 
models. We do not aim for an accurate phenomenology of the nucleon and the A 
resonance. It is known from the experience gained in local models that, in addi- 
tion to the mean-field approximation applied in our work, many other effects can 
and should be included: the projection of the center-of-mass motion and angular- 



momentum, ||48|j49| , |50| , |51| , |52| , the rotational 1/Nc corrections ||53|j5^ , inclusion of 



other degrees of freedom such as vector mesons [Pl,pH] or color-dielectric fields 



3^. These may considerably alter and improve the results. We stress that the 
methods which we describe in this paper are applicable to all models with separa- 
ble nonlocal four-fermion interactions provided that it is possible to perform an 
analytic continuation in the four- momentum close to zero, as explained in Sect. 

w. 



The existence of stable solitons in nonlocal models has been briefiy reported in 
Refs. [^,^. In the present work we describe in greater detail how the soliton is 



constructed and how the valence orbit is defined (Sect. |2.4| ). The determination 



of the model parameters is described in Sect. ^ The physical properties of the 
soliton and its energetic stability are analyzed in Sect. ^ An important and novel 
feature of this paper concerns the gauging of the nonlocal model and the con- 
struction of Noether currents. In Sect. |^, we apply the technique of path-ordered 
P-exponents and we obtain general expressions for currents in the soliton back- 
ground. Interesting results follow: the properties of the soliton, which are related 
to currents at zero momenta (charges, qa-, and in general ra-point Green functions 
with vanishing momenta on external legs) do not depend on the choice of the 
path in the P-exponent, and are thus universal. On the other hand, observables 
which probe the soliton at non-zero momenta (form-factors, magnetic moments) 
do depend on the gauging prescription. However, in the cases we have explored, 
the contributions of the nonlocal parts of currents are small, in particular for 
sufficiently large solitons. In Sect. || we apply two prescriptions to compute ob- 
servables: straight-line paths in the P-exponents and the weak-nonlocality limit, 
where the energy scale of the regulator is assumed to be parametrically much 
larger than the other scales in the problem. 

The many-body techniques applied in our work are inherited from the experience 
acquired with local models. In particular, the nucleon is derived from the hedge- 
hog soliton by cranking. We calculate the electric isoscalar radius, the isovector 
magnetic moment an the isovector magnetic radius. These quantities do not de- 
pend dynamically on cranking |^^. They are much simpler to evaluate numeri- 



cally, involving single spectral sums, than the moment of inertia, isovector electric 
radius, or the isoscalar magnetic properties, which involve double spectral sums. 

Appendix ^ contains a detailed account of the gauging method. Appendix ^ 
describes the construction of general Noether currents. Explicit forms for the 
straight-line P-exponents are given in App. and the forms obtained from the 
weak-nonlocality limit are given in App. ^ Expressions needed for the evaluation 
of observables are displayed in App. 0. 



2 The model 



This section describes the model used to calculate the solitons. We derive the 
basic formulas, such as the Euler-Lagrange equations, and the expressions for 
the baryon number and energy. We discuss the construction of the valence orbit. 



which is non-trivial in nonlocal models 37 . 



2.1 The action 

The model is defined by the Euchdean action 
I = j d^x q {-id^-ff, + m)q- — {^ipTaip 



(1) 



where q (x) is the quark field, m the current quark mass, 7/^ = 7'^ = {iP, 7) are 
the antihermitian (Euclidean) Dirac matrices, and Fq = 1, F^ = i^^Ta define the 
coupling in the scalar and pseudoscalar qq channels. We work with Nf = 2 fiavors. 
The coupling constant G has the dimension of inverse energy. The delocalized 
quark field, ip (x), is related to the quark field q{x) by a regulator r, which is 
diagonal in the momentum space: 

ij (x) = {x\ij) = {x \r\ q) = -^Yl rke'''"'qk , (2) 

k 

where f2 = / d4X is the Euclidean space-time volume. 



Our calculations use either a Gaussian (as first considered in Refs. [^,^) or a 
monopole [^] form for the regulator: 

_ fe2 1 

Tk = e 2A^ (Gaussian) , r^ = r^- (monopole) . (3) 

The action (|l]) is easily cast into the form 
I = -F:J2^k {kf,-ff, + m) qk 

k 

~^ Z! ^k^+k^M+kA rk^rk^rk^rk^ {qkj^aqk2) {qkj^aqk^) ■ (4) 

kik2k:jk4 

The regulators appearing in the four-quark interaction have a separable form in 
momentum space. The separability of the quark interaction is also present in the 
instanton-liquid model 0,|^. 

Soliton calculations are more easily performed by using the equivalent bosonized 
form of the action (|l|): 

I {S,P) = -TT\n{-id,-f, + m + r<!>r) + ^ J d,x [S^ + P!) , (5) 

where 

^{x) = Six)+l^,Paix)Ta (6) 



is the local cliiral field which is the dynamical variable of the system. The trace 
in Eq. (^ is over color, flavor, Dirac indices, and space-time. The chiral fields 
represent the following expectation values of bilinear forms of the quark fields: 

S (x) = —G'^{iI) (x) il> {x)) = —G^ {q \r\ x) {x \r\ q) , 
Pa {x) = -G'^{i) (x) i-^^Tai) (x)) = -G^ (g \r\ x) i'ji^Ta (x \r\q) . (7) 



Since stationary solitons involve time-independent S and P fields, it is more useful 
to express the action in terms of the Dirac Hamiltonian 

a ■ V 
h = — : h /3m + /3r $ r, (8) 



bearing in mind, however, that the regulator makes the Dirac Hamiltonian de- 
pendent (although diagonal) on the energy variable. The action (^ becomes 

/ = -Tr ln(9. + h) + ^Jd,x {S' + Pf) , (9) 

where r is the Euclidean time variable. 



2.2 The hedgehog soliton 



Our use of the model action (|^) is the same as in other hedgehog soliton cal- 
culations: we treat the S and P fields classically, thus keeping the leading-order 
contribution in the number of colors, Nc- Interestingly, in nonlocal models the 
next-to-leading-order corrections are found to be surprisingly small in the vacuum 
sector |2|j5|j6gj61|j6gj6|. 



We seek a stationary point of the action, or, more precisely, a minimum of the 
energy for time-independent chiral fields with a hedgehog shape 

S{X)=S{\X\), Pa{x)=XaP{\x\), (10) 

where Xa = Xa/|x|. Asymptotically, far from the soliton, the fields recover their 
vacuum values S (x) = M and Pa (x) = 0. We shall refer to the vacuum value M 
of the scalar field as the constituent quark mass (see the discussion of what we 
mean by the quark mass in Sect. jOP . The soliton represents a bound state of 
Nc = 3 quarks occupying a valence orbit with the grand spin and parity G^ = 0"^, 
where the grand spin G = J + / is the sum of the total angular momentum and 



isospin p0| , |32| , together with a Dirac sea which is polarized by the hedgehog field. 



2.3 The energy- dependent quark orbits 

The regulator acts as a differential operator, r = r ( —d^ — V^ ) . It introduces a 
nonlocal interaction between the quarks and the chiral fields. For time-independent 
chiral fields $, defined in Eq. (|^), the Dirac operator is diagonal in the energy 
representation, 

(dr + h) \uj) = (iu + hito^)) \uj) , (11) 

where h (lj^) is the energy-dependent Dirac Haniiltonian: 

h{uj^) = 1^^ + (3m + f3r{u^ - V') $ r{uj^ - V'). (12) 

i 

The Kahana-Ripka basis [^] is convenient to diagonalize the Dirac Hamilto- 
nian h, since the regulator r is diagonal in this basis. Indeed, using the notation 
\skljGG3) for the Kahana-Ripka basis states, we have r{uj'^ — V^) \skljGG3) = 
r{uj'^ + k"^) \skljGG3). In this notation, k is the radial momentum of a quark 
quantized in a spherical box, / is the orbital momentum of the upper component, 
j is the total angular momentum of the quark, and G is the grand spin, obtained 
by adding j and the isospin 1/2. 



We diagonalize the Dirac Hamiltonian ( p!2|) for each value of cj, thereby obtaining 
the energy-dependent quark orbits: 

/i(cu2)|A^) = eAp)|A^). (13) 

Observables can then be calculated in terms of the quark orbits \Xuj). A technical 
complication, compared to earlier calculations which used local models, is the 
presence of an additional integral over u in the expressions for observables, which 
has to be carried numerically. 



2.4 The energy, baryon number, and valence orbit of the soliton 

We can use the basis \lj) |A^) to express the action (|^) in the form: 

J = - ^ In (.^ + e, (^2)) + ^Jd,x {S' + Pj) . (14) 



The sum over A^ includes color. In the zero-temperature limit and for time- 
independent fields S and P the energy of the system is equal to 



^ = - / ^ E In 0^ + eA(c.^)) + ^Jdsx {S' + Pf ) - E^^ , (15) 



where -Evac is the vacuum subtraction, obtained by setting S (x) = M and 
Pa (x) = 0. This subtraction ensures (together with the regulator) that the energy 
remains finite. We integrate the first term by parts. The boundary term cancels 
out when the vacuum energy is subtracted, and we find 



oodoo y-^ 1 — i 



dex 



^=r-^Y.:—^)^^Ji^-(s'^pi)-E,,.. (16) 



The determination of the baryon number of the soliton is not obvious because 
the regulators appearing in the action (|l]) prevent us from performing a canonical 
quantization of the quark field. In spite of this, the baryon number turns out to be 
correctly quantized. A simple way to see this consists in constructing the Noether 
current associated to the abelian gauge transformation q (x) — > e~'^^^'^'^q (x) of the 
quark fields |jl3[ . For the sake of calculating the baryon number we assume that 



?7 (r) depends only on time (complete Noether currents are derived in App. P). 
The quark loop term of the action is transformed to 

Jg (r/) -^ -Tr In e'"^ {dr + h) e''" = -Tr In (-ir/' + e'^^he''") , (17) 

where rj' = drj/dr. In the Noether construction of the baryon number extra terms 
arise due to the nonlocal regulator r, which does not commute with rj. The Dirac 
Hamiltonian is diagonal in uj and rj is diagonal in r. To first order in rj' we can 
use the commutator expansion 

e^^/ie"*" = h + i[r],h] = h- 7]'h', (18) 

where h' = dh/duj. Therefore 

/, (r^) = /, (r^ = 0) - Tr In [d^ + h - iri (1 - ih!)) , (19) 

and the baryon number is 

where the factor XjNc is due to the fact that the quark carries baryon num- 
ber 1/iVc. Using the spectral basis \ijS) \\J) to evaluate the trace, we obtain the 



following expression for the baryon number of the system: 

B = - /^V \ t . - (21) 



The presence of the factors 1 — ^-Jj is crucial for the quantization of the baryon 
number. They appear in the calculation of all observables related to Noether 
currents, such as the baryon density, qa, and magnetic moments, as shown in 



■ de 
du 

propagator 1/ {dr + h). Indeed, in the vicinity of a pole at ujq we find 
1 1 1 



App. ^. The factors 1 — i-^ are the inverse residues of the poles of the quark 



■ (22) 



uj-iex{uj) uj-uo-i^^ 



UJO 



{UJ - UJo) ^ - ^0 ( I _ I '^eA('^) \ 



The position of these poles corresponds to the "on-shell" single-quark energies. If 
we were able to deform the integration contour of the energy variable uj such as 
to close it at infinity above or below the real axis, the energy would become the 
sum of the on-shell quark energies, and each pole would contribute a factor of 1 
to the baryon number (strictly, a factor l/Nc-, but we construct colorless solitons 
by placing Nc quarks into each orbit). The baryon number is thus effectively 
quantized despite the fact that we cannot perform a canonical quantization of the 
quark fields. 

We now construct the valence orbit in the nonlocal model. In our calculations the 
spectrum of the Dirac orbits lA^^) in the hedgehog field displays a similar pattern 
for all values of the Euclidean energy variable uj. It consists of a Dirac sea, which 
is the set of negative energy orbits and a set of positive energy orbits, which are 
separated from the Dirac sea orbits by a well defined energy gap. Within this 
gap, there appears a bound positive energy (/ + J)^ = 0"*" orbit, which we call 
the valence orbit and which is well separated from both the positive energy and 
negative energy Dirac sea orbits. We do not observe crossings between these sets. 



As a result, the baryon number (|2l|) of the polarized Dirac sea is the same as the 
baryon number of the vacuum, namely 0. In order to obtain a soliton with the 
baryon number equal to 1, it is necessary to deform the integration contour in 
Eq. (|21| ) as displayed in Fig. 1, such as to enclose the valence pole. The valence 
pole is the solution of the equation 

u = iey{u;'^). (23) 

The energy e^ (cj^) is the eigenvalue of the Dirac Hamiltonian h (cj^) correspond- 
ing to the valence orbit. Equation (^) is a non-linear equation for uj and the 
solution is written as icvai- It can be found numerically without difficulty or am- 
biguity. We refer to Cvai as to the valence energy, and the corresponding valence 



& 



Fig. 1. The contour in the integral over the Euchdean energy variable, uj, encircling a 
positive-energy valence orbit. 

state, |val), satisfies the equation h (— e^^^i) |val) = evai|val). When the integration 
path of u is deformed so as to include the valence orbit, as shown in Fig. 1, the 
soliton acquires baryon number B = 1, and its energy is equal to 



E = iV.evai + 



ujduj 
~2^ 



■, dex 
duj 






^Jd,x {S' + P',)-E^ 



(24) 



The first term is the valence orbit contribution, the second term, with the in- 
tegration carried along the real axis in the u complex plane, is the Dirac-sea 
contribution. The third term is the energy of the meson fields. 

In addition to the valence pole, located on the imaginary Euclidean energy axis 
in the vicinity of 0, there are many other poles in the complex-energy plane. This 
feature is present in nonlocal models, and also in certain local models, e.g. those 
with the proper-time regularization |]6^. Fortunately, for the chosen regulators 
@ the other poles lie far away from the origin on the scale A, hence they do not 
interfere with the valence orbit, and the deformation of the cu-integration contour 
is well-defined. 



The pole defined by Eq. (^) occurs at an imaginary u and therefore at a neg- 
ative u"^. The determination of the valence orbit requires therefore an analytic 
continuation of the regulator to negative arguments because —e^^^ + k'^ can be- 
come negative for low-enough fc^. Such analytic continuation causes no problem 
for the forms (H) of the regulator, but it prevents us from using the regulator 



derived from the instanton model of the vacuum 14 , because the latter has a 



branching point at fc^ = 0. The analytic continuation is only required for small 
negative values of k"^ = uj"^ + k"^, considerably smaller than the nonlocality scale 
A^. Indeed, Cvai lies within ±M, and, on the average, k"^ = —el + k"^ remains very 
small relative to A^. 
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2. 5 The Euler- Lagrange equations for the soliton 



The static meson fields are determined self-consistently by solving the Euler- 
Lagrange equations derived by minimizing the energy (|T^) with respect to vari- 
ations of the fields S and P\^ 



oo 



A^cZvai(valr x)/5 xrval)+ / 7^2^ ■ — ■ ^^r — , (25) 

-00 -^^ 

A^cZvai(val|r|x)?/375ra(x|r|val) + / 7— V" 

J Ztt , 



——- = A^cZvai val r X ^/375r„ X r val + ttI^ ■ — ^ Tr. 



where |val) is the valence orbit, and the residue factor is 

-1 
, (ievai(^) 

Zval = 1 - « 



duo 



(26) 



<^=«eval/ 



Because the energies e\ {u"^) and the wave functions \\J) depend only on cj^, we 
may symmetrize the integrands with respect to uj and rewrite the Dirac-sea parts 
of Eqs. (^5|) in the manifestly real form 



-00 
00 



du ^ ex (g;^) {\u,\r\x)[3{x\r\\J) 

2vr t ^^ + e! (^^) ' ^ ^ 

27r f- (x;2 + ei (u;2) 

-00 -^IJ A \ / 

The soliton is calculated iteratively. An initial guess is made for the fields S and 
P. In terms of these, the quark orbits (|TBp are calculated for different values of uj. 
The fields S and P are then recalculated using (^), and so on until convergence 
is reached. 



3 The vacuum sector 



The vacuum sector of our model, analyzed in Refs. pO| , piiP2| , is used to constrain 
the model parameters, namely the strength of the quark-quark interaction G, 
the cut-off A, and the current quark mass m. We fit the pion decay constant. 



^ The Euler-Lagrange equations in Ref. |3^ had a typographical error by which the 
residue factors Zvai appearing in Eqs. ([25|) were omitted. 
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fn = 93 MeV, and the pion mass, m.,r = 139 MeV. This leaves one undetermined 
parameter, which we choose to be M , the vacuum value of the scalar field 5*. 
Admittedly, it is somewhat abusive to claim that the model depends only on 
three parameters. It is only true to the extent that the regulator r^ is expressed 
in terms of one parameter, such as the cut-off appearing in (^. The model depends 
in fact on the shape of the regulator and furthermore, as we have seen in Sect. 
27^, the construction of the valence orbit relies on the possibility of the analytic 



continuation of the regulator to (small) negative momenta k . 

For smooth regulators it is found that quantities such as masses, /vr, or 771^^, which 
would diverge logarithmically with the cut-off, do not depend very much on its 
shape. When one of these quantities is fixed, /^r for example, the masses of mesons 
and of solitons do not depend very much on the form of the regulator. However, 
quantities, such as the quark condensate, which diverge quadratically with the 
cut-off, are more sensitive to the shape of r^ and different values are obtained with 
different shapes even when f^^ is kept fixed [^] . This can be seen from Table |I[ In 



the whole range of values 300 < M < 600, the cut-off of the Gaussian regulator 
remains 1.06 times larger than the cut-off of the monopole regulator, within 1% 
accuracy. The two regulators have the same low-fc behavior, therefore the shape 
dependence is roughly a 6% effect on quantities such as f-,, or m.„. 

One may rephrase the above statements as follows: for any given M we fit A for 
each regulator (Gaussian and monopole) in such a way that /^ is fixed to 93 MeV. 
Then, the regulators (^, each with its own A, are very similar functions of the 
k"^ variable up to A;^ ~ 7 fm~ when M = 300 MeV, and k"^ :i^ 3.5 fm~ when 
M = 600 MeV. Therefore, the quantities dominated by low values of k'^ depend 
only weakly on the regulator. 



3. 1 The quark propagator in the vacuum 



In the vacuum sector, the fields acquire the values S = M and P = 0, and 
the inverse quark propagator is diagonal in momentum space. In the chiral limit 
m — > 0, it is equal to k^'y^ + rlM. Poles of the quark propagator occur when 
k'^ = — M^r'*(/i;^), where k'^ = u'^ + k"^. The solution of this equation can be 
written in the form k"^ = — M^ , where Mq satisfies the equation: 

M^ = M'^r^ {-Mg) . (28) 

This non-linear equation for M^ has, in general, solutions which are scattered in 
the complex k"^ (or M^) plane. If M/A is small enough, a solution occurs on the 
real axis of the M^ plane and such a pole represents an on-shell free quark with 



12 



Regulator 


M 


A 


m 


(k-^)'/' 


1/G 


(^G^^GD^/^ 




[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 




300 


760 


7.62 


-222 


182 


327 




350 


627 


10.4 


-200 


140 


310 


Gaussian 


400 


543 


13.2 


-185 


113 


298 




450 


484 


15.9 


-174 


94 


287 




600 


380 


24.2 


-151 


61 


268 




300 


718 


3.96 


-276 


204 


347 




350 


590 


5.24 


-252 


159 


334 


monopole 


400 


509 


6.44 


-235 


130 


323 




450 


452 


7.56 


-223 


110 


315 




600 


352 


10.5 


-200 


74 


294 



Table 1 

The vacuum properties for the two regulators considered, and for various values of the 
constituent quark mass M. For each case the values of F„ and m-,r have been fixed to 
their physical values. 



a mass equal to Mq, which we call the vacuum on-shell quark mass. This on-shell 
quark mass Mq can be considerably larger than M, which we call the constituent 
quark mass. As we shall see in Sect. §, it is the on-shell quark mass Mg, and not 
the constituent quark mass M , which determines the stability of the soliton (see 
the discussion of Fig. 2). When the energy (|2^) is greater than N^Mg, the soliton 



is not formed, which means that the Euler-Lagrange equations (p5|) do not have 
a localized stationary solution. The soliton can thus be viewed as a bound state. 
The stability of the soliton is discussed at the end of Sect. H. 



3.2 The gap equation and the condensates 



When discussing the model in the vacuum sector, it is much simpler to use the 
Lorentz-invariant form (^ of the action. In the vacuum sector, a translationally 
invariant stationary point of the action exists with S = M, P^ = 0, where M is 
the solution of the equation: 



^ = 4iVciV^ 



ddk 



{2ny P + (m + rlMY 



(29) 



13 



Equation (p9D is traditionally called the "gap equation" . It is the Euler-Lagrange 
equation expressed for a translationally-invariant system (without valence quarks). 
For a given current quark mass m (which is determined by fitting m^), the gap 
equation (^9]) relates M to the interaction strength G. In this work we use it in 
order to eliminate the parameter G in favor of M. 



The quark condensate (qq) = (uu + da) can be obtained from the action (§). In 
the chiral limit it is equal to 



I dl 



ANcNf 



m=0 



ddk 



(27r)^ P + (rlM) 



2- 



(30) 



In the instanton model of the QCD vacuum the gluon condensate can be expressed 
in terms of the constituent quark mass and the four-fermion coupling constant 



a. 



IT 






32Nr 



ddk 



rJM^ 



8M^ 



(27r)4 p + (riM)' NjG^' 



(31) 



The numerical values are listed in Table |I|. The estimate for the gluon condensate 



inferred from the QCD sum rules |67,68[ is (—G 



360 ± 20 MeV, while Ref. 
169|] gives the value 386±10 MeV. These estimates, when compared to the numbers 



Al/4 



in Table |I|, favor lower values of M in our model. 



3.3 The pion mass and the pion decay constant 



The inverse pion propagator Kp^{x,y) in the vacuum can be deduced from the 
action (^): 



xa 



Kp' 



yb 



5^1 {S, P) 
5Pa{x)5P,{y) ■ 



(32) 



Explicit expressions for the meson propagators, /^r, and 771,^, calculated with the 
present model, can be found in Refs. [p^pO| , pT| , p^ ,|6^ and we shall not repro- 
duce them here. We only specify the expressions used to determine the model 
parameters. The pion propagator is diagonal in momentum and flavor space: 
q'h) = SabSqq'Kp^ {(f) ■ To lowest order in q^ and m it acquires the form 



qa 



Kp' 



Kp' {q') = Z^ {q' + ml) , 
where the inverse residue of the pion pole is 

Z^ = 2N,Nf f (M) 



(33) 



(34) 
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and where the function / (M ) is 



/(M) 



»,4 _ 1.2^2 drl , ,4 ( drl 

d.k ^k f^ r^^ + ft — 



kdk2 



dk2 



(27r) 



{k^ + riM-^y 



The pion decay constant, fj^, is then given by the expression |l^ , pO| ,pT 
/2 = z^M' = 2N,NfM^ f (M) , 



and the pion mass equals to ||l^ , pO| , pT 

n 2m 



m^ 



M f{M) 



g'{M), g'{M) 



dik 



(27r)^ fc2 + rtM^ ' 



(35) 



(36) 



(37) 



From Eq. ( pO]) we can see that, in the chiral hmit, the quark condensate is (qq) 
-ANcNfMg' (M). Using (||), Eq. (|^ can be cast into the form 



m. 



m {qq) 



(38) 



which is the Gell-Mann-Oakes-Renner relation, requested by the constraints of 
the chiral symmetry. 



The model parameters M , A, m, and G are listed in Table [I|. We notice that 
when M increases, the cut-off A decreases. This occurs because the pion decay 
constant is kept fixed at the value /^r = 93 MeV. For large A, the quark condensate 
increases linearly with M and decreases quadratically with A. The net result is 
a slow decreases with M. The coupling constant G of the attractive quark-quark 
interaction increases with M. This is why solitons are bound only when M is 
large enough, as discussed in Sect. §. At large values, M > 450 MeV, the cut-off 
becomes embarrassingly small as compared to the other scales in the problem. The 
error of taking the leading expressions in m extends from 2 % at M = 300 MeV 
to 10 % at M = 600 MeV. This estimate is based on the difference of the integrals 
(|30D , (^1]) , (p5D and (|3^) in the chiral limit with m from Table |I|. 



4 Properties of the soliton 



The soliton is calculated self-consistently by solving iteratively the Euler-Lagrange 
equations ( ^3]) for the meson fields and the Dirac equation for the quark orbits 
(|13D, as described in Sect. ||. The convergence is fast except for very low values of 
M. Because of the presence of the regulator which cuts off very high momenta, it 
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Fig. 2. The energy of the sohton (bold sohd hne), Nc times the free-space on-shell quark 
mass (thin sohd hne) , and the valence contribution to the soliton energy (dashed line) , 
plotted as functions of the parameter M. The dotted line shows the cut-off parameter, 
A, fitted to obtain /^r = 93 MeV. All quantities in MeV. The Gaussian regulator is 
used. 

is not necessary to use a very large basis as in similar calculations in local models. 
We have performed the calculation with two shapes (|) of the regulator, Gaus- 
sian and the monopole (^. The soliton properties are sensitive to low values of 
k"^, such that, according to the discussion in the beginning of Sect. |^, the results 
depend very weakly on the shape of the regulator [^ (see Table ^) [^. 

Figure 2 shows the energy of the soliton obtained with the Gaussian regulator. 
All quantities are plotted as a function of the constituent quark mass, M. The 
energy of the soliton is a slightly increasing function of M while the energy of the 
valence orbit flattens out. The soliton energy varies from about 1100 to 1250 MeV 
when M increases from 300 MeV to 450 MeV. These seem to be reasonable values 
for a soliton which is to describe the nucleon, because the energy is expected to 
decrease by about 200 MeV when the center-of-mass energy is subtracted in a 
suitable projection method |4^j49|j50|j5T|j52 1 . 



^ This is true for regulators which, as a function of fe^, have non-zero slope at the 
origin. No stable solutions are found for regulators for which this derivative is zero. 
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Fig. 3. Self consistently determined fields for various values of M, plotted as functions 
of the radial coordinate x. The Gaussian regulator is used. 



The curve labelled 3Mg gives Nf. = 'i times the value of the on-shell quark mass 
Mg in the vacuum. It is a solution of Eq. (p8|). The curve terminates when M ~ 

309 MeV, or, more precisely, when ^ ~ J-^ = 0.43. (The exact value depends 
on the shape of the regulator and on m.) Indeed, for larger values of ^, Eq. 
(pSD no longer has real solutions and the "on-shell" quark mass wanders off into 
the complex plane. This has been claimed to be related to quark confinement 
[|6l,|l8| , although the nature of this relation is far from clear. The on shell valence 
orbit in the background hedgehog field, defined by Eq. (^), exists even when the 
model parameters prevent the occurrence of an on-shell quark pole in the vacuum 
background field. 

A bound state of quarks occurs when the soliton energy E is lower than the 
energy of A'^c = 3 on-shell quarks, that is, when E < NcMq. Figure 2 shows that 
a bound state of quarks only occurs when M exceeds a critical value of 276 MeV 
{i.e. when the coupling constant G > 4.12/A ). Beyond the critical value of M 
where Mg becomes complex, a stable solution continues to existQ. A very similar 
behavior is found for the monopole regulator. 

Figure 3 shows the radial shapes of the hedgehog chiral fields S (x) and P (x) 
in units of M. They are the solutions of the Euler-Lagrange equations (pSj). We 



^ Since complex poles cannot be asymptotic states of the theory, the soliton (or any 
other object) can never decay into such states. 
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M = 600 MeV 
M = 450 MeV 
M = 350 MeV 
M = 300 MeV 



_] I I L 



12 3 4 

X [fm] 

Fig. 4. Effective squared quark mass in tlie soliton plotted as a function of the radial 
coordinate x. The Gaussian regulator is used. 

can see that, in the reasonable range 350 MeV < M < 450 MeV, the chiral 
field deviates significantly from the chiral circle. Only at excessively high values, 
M > 600 MeV, does the chiral field remain close to the chiral circle. This is a 
new dynamical result. In previous soliton calculations, which used the Nambu 
Jona-Lasinio model with proper-time or Pauli-Villars regularizations [0,^, the 
soliton was found to be unstable (the energy was unbounded from below) unless 
the fields were artificially constrained to the chiral circle |HO|,ffT|. Our model is free 



from this restriction. In our soliton, the self-consistent pion field is considerably 
smaller, as compared to previous calculations. In a sense, it is midway between 
the Skyrmion (where the chiral field is on the chiral circle) and the Friedberg-Lee 
soliton ||7I| , |72| (where the pion field is zero). Deviations from the chiral circle are 
further illustrated in Fig. 4 which shows the values of S*^ (x) + P^ (x) in units of 
M^. The curve would remain constant and equal to 1 if the fields remained on 
the chiral circle. Another way to phrase the behavior displayed in Fig. 4 is to say 
that the chiral symmetry is partially restored in the center of the soliton. 

Figure 5 shows the upper (m) and lower (t>) quark components for the valence 
orbit for various values of M. We note that the solution shrinks as M is increased, 
however, beyond M ~ 450 MeV the effect saturates. 

Figure 6 illustrates the stability of the soliton, composed of 3 quarks, with respect 
to its fragmentation into solitons formed with 1 or 2 quarks. Due to the lack of 
confinement, such solutions formally exist in the model. The solid line, labelled 
\Esoi (3^), shows the soliton energy per quark. The dashed and dot-dashed lines 




M = 600 MeV 
M = 450 MeV 
M = 350 MeV 
M = 300 MeV 



Fig. 5. The upper (u) and lower (v) quark components for the valence orbit for various 
values of M, plotted as functions of the radial coordinate x. The Gaussian regulator is 
used. 

show the energy of the single-quark and two-quark soliton, respectively. The curve 
labelled Mg is the on-shell mass of the quark in the vacuum. We conclude from 
Fig. 6 that the 3-quark soliton is stable against the breakup into solitons with 
a lower number of quarks. Similar results have been found in the linear sigma 
model with valence quarks ||73|. Note that the Pauli principle prevents placing 
more than Nc = 3 quarks into the 0"*" valence orbit. 



5 Noether currents 



The Noether currents in nonlocal models acquire extra contributions due to the 
momentum-dependent regulator. As mentioned in the introduction, the trans- 
verse parts of Noether currents are not fixed by current conservation and their 
choice is not unique. Any prescription becomes an element of the model building. 
The problem of ambiguous transverse currents has been known for a long time in 



connection with meson-exchange currents pO| , ^ . An elegant way of gauging the 
nonlocal model is to use path-ordered P-exponents, and we choose this technique 
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Fig. 6. The energy per quark for the sohton with three valence quarks (bold solid line), 
one valence quark (dashed line), and two valence quarks (dashed-dotted line), plotted 
as functions of the parameter M. The thin solid line shows free-space on-shell quark 
mass Mq. All quantities in MeV. The Gaussian regulator is used. 

to construct the Noether currents. 

The details are shown in Appendices 0-^ In this section we make some general 
remarks: 

- The ambiguity in the choice of Noether currents is attributed to the freedom 
of choosing the path in the P-exponent. 

- The Noether currents associated to symmetries are conserved. 

- The properties of solitons involving zero-momentum probes, such as the baryon 
number and qa, do not depend on the chosen path in the P-exponent. They 
are thus defined unambiguously. 

- Soliton radii, magnetic moments, form factors, do depend on the path, and 
hence they are not uniquely determined. We find, however, that the path- 
dependence is weak in the weak-nonlocality limit, i.e. in the case where the 
soliton scales are much smaller than the nonlocality scale, A. This is the case 
for sufficiently large solitons. 

Using a straight-line path in the P-exponent, we derive in App. |C| the following 
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Regulator 


M 


fival 


-C/sea 


-C/mes 


^sol 




[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 




300 


298 


2420 


-2227 


1088 




350 


285 


1773 


-1450 


1180 


Gaussian 


400 


279 


1494 


-1102 


1228 




450 


275 


1339 


-905 


1260 




600 


272 


1117 


-619 


1313 




300 


289 


3008 


-2790 


1084 




350 


275 


2201 


-1850 


1176 


monopole 


400 


266 


1835 


-1407 


1227 




450 


260 


1628 


-1147 


1261 




600 


248 


1332 


-753 


1321 



Table 2 

Contributions to the soliton energy calculated with various regulators. 

form for the nonlocal contributions to Noetlier currents: 

duj 



■NL.straight ^ g\ 



/ —J2h^3xd3yda6{z-x-a{y-x)) x 

A^ 



(A^|f) {x\r^{uj) \y) {y\\aP^r{u;)\X^) 
iuj + ex iuj) 



+ h.c, 



(39) 



where A^ stands for 1/iVc, ra/2, or ^^Tajl in the case of the baryon, isospin, or 
axial currents, respectively. The parameter a describes the straight-line integra- 
tion path. The space integrations reflect the nonlocality. The form (|39D looks 
somewhat complicated. However, moments of currents i^gA-, radii, or magnetic 
moments) are very easily evaluated because, for those cases, the a integration 
and one space integration can be performed trivially (see App. |^ for examples) 
leaving only one space integral, as in the local contributions. 

A simpler expression for the Noether currents can be derived in the weak-nonlocality 
limit [0. As shown in App. 0, we have: 



j^'L"'"^^' (^1 



TT-i: — — + h.c. 



277 



tuj + ex [uj) 



(40) 



In the calculation of observables, we single out the valence and Dirac-sea parts 
in Eqs. (^,^) in the usual way, as shown in App. ^ 
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6 Observables 



As is well known, the nucleon quantum numbers (spin, flavor) need to be pro- 
jected out of the hedgehog soliton in order to calculate observables. In the large- iVc 
limit this can be achieved by cranking []57| , [71| . The observables fall into two cat- 
egories, according to whether they are dynamically-independent or dynamically- 
dependent on cranking. The former ones include qa, isoscalar-electric, and isovec- 
tor magnetic properties, while the latter ones include isovector electric and isoscalar 
magnetic properties. Observables which are dynamically-independent of cranking 
are simpler to evaluate, since they involve only the expectation value in the hedge- 
hog state. Technically, they lead to single spectral sums over the quark orbits. 
Quantities which are dynamically-dependent on cranking involve double spectral 
sums, and they are difficult to evaluate even in local models |]l2i|39| . The presence 



of nonlocality adds extra difficulties: an integration over the energy variable, and 
the appearance of nonlocal contributions. Because of these difficulties, and since 
the aim in this paper is primarily to develop a consistent approach to calculate 
observables in the presence of nonlocal regulators, we restrict our calculations to 
the cranking-independent observables. 



6.1 The isoscalar electric radius 



The results for the mean squared isoscalar electric radius, calculated with the 
help of the formulas given in App. [E.l| , are displayed in Table |]for various values 
of M. The local contribution and different nonlocal contributions (see App. [E.ll 
for their meaning) labeled NL(A), NL(B) and NL(C), are listed separately for 
the valence orbit and for the Dirac sea. As discussed in the previous section, the 
radius depends on the chosen path in the P-exponent, and the result is not unique. 
We give the results for a straight line prescription and for the weak-nonlocality 
approximation . 

The expressions for the nonlocal terms contain a derivative of the regulator which 
produces a factor 1/(2A^), which suppresses the nonlocal contribution, as com- 
pared to the local one. Since A decreases with increasing M , the nonlocal terms 
become more and more important for larger M and the difference between the 
two prescriptions of evaluating the Noether currents increases. The difference is 
reasonably small for the physically relevant values of M. The soliton is weakly 
bound for the values of M just above the threshold and hence very large. Its size 
decreases and reaches the minimum around M = 450 MeV. Above this value, the 
radius starts increasing again. This is due to the fact that A is very small and 
the size becomes proportional to the inverse A. 
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Fig. 7. Various contributions to the baryon density (multiplied by Attx'^) for 
M = 350 MeV and the Gaussian regulator, plotted as functions of the radial coor- 
dinate x: the local (solid bold line) and nonlocal (solid thin line) valence contributions, 
and the local (bold dashed line) and nonlocal (dashed thin curve) sea contributions. 
The sea contributions are multiplied by a factor of 100. 



The nonlocal contribution from the valence orbit considerably reduces the isoscalar 
electric radius, up to 20% for low values of M. This is consistent with the fact 
that the inverse residue, Zvai, and hence the integrated local density, is greater 
than 1. The corresponding nonlocal density is negative for all \x\ (see Fig. 7). 
We can also see that the weak-nonlocality limit works better for lower values 
of M. This is natural, because the weak-nonlocality limit may be viewed as the 
leading-order term in an expansion in powers of l/(i?A), where R describes the 
soliton size. This is why the B terms, carrying more derivatives of r, are much 
smaller than the A and C terms. By the same argument, in the weak-nonlocality 
limit A = C. We note finally that our numbers for (r^)baryon are larger than the 
experimental value of 0.62 fm^. 

Different contributions to the baryon (isoscalar electric) density are displayed 
in Fig. 7. As a numerical check we have verified that the sum of the valence 
contributions integrates to 1, while the sum of the sea contributions integrates to 
0. 
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M [MeV] 


( 


2\ 
^ / baryon 


300 


350 


400 


450 


600 


val 


L 


2.209 


1.270 


1.057 


0.991 


1.062 




NL (A) 


-0.446 


-0.228 


-0.188 


-0.192 


-0.328 




NL (B) 


-0.039 


-0.039 


-0.040 


-0.041 


-0.061 




NL (C) 


-0.435 


-0.194 


-0.119 


-0.074 


0.021 




straight line 


1.761 


1.052 


0.897 


0.852 


0.898 




weak NL 


1.773 


1.075 


0.937 


0.917 


1.083 


sea 


L 


0.0050 


0.0070 


0.0080 


0.0080 


0.0075 




NL (A) 


0.0005 


0.0012 


0.0020 


0.0022 


0.0030 




NL (B) 


0.0001 


0.0001 


0.0002 


0.0002 


0.0003 




NL (C) 


0.0003 


0.0010 


0.0015 


0.0018 


0.0026 


total 


straight line 


1.766 


1.060 


0.907 


0.862 


0.909 




weak NL 


1.778 


1.083 


0.947 


0.927 


1.094 



Table 3 

Various contributions to the isoscalar electric mean square radius. L denotes the local 
contribution and NL(A,B,C) different nonlocal contributions defined in App. E.l. The 
total result is displayed for the straight path and in the weak-nonlocality approximation. 
The Gaussian regulator is used. 



6.2 QA 



The results for qa, evaluated with the expressions given in App. [E.2| , are displayed 
in Table ^. The sea contribution remains small for all values of M. The nonlocal 
terms increase with increasing M, as expected, yielding, together with the local 
piece, almost a constant value of qa over a wide range of M. The values, ranging 
between 1.1 and 1.15, are somewhat smaller than the experimental value of 1.26. 



6.3 Isovector magnetic moments and radii 



The results for the isovector magnetic moment, evaluated with the expressions 
given in App. |D.3| , are shown in Table |^. As in the case of qa, the nonlocal terms 
increase with increasing M, and the total value is almost constant over a wide 
range of M. The values are lower than the experimental value 4.71. 
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M [MeV] 


9A 


300 


350 


400 


450 


600 


val, L 

val, NL 


L047 
0.022 


0.922 
0.069 


0.861 
0.119 


0.819 
0.170 


0.737 
0.309 


sea, L 
sea, NL 


0.050 
0.032 


0.067 
0.053 


0.064 
0.062 


0.058 
0.065 


0.037 
0.062 


total 


1.151 


1.112 


1.106 


1.112 


1.146 



Table 4 

Different contributions to qa calculated from the current. L and NL denote the local 

and nonlocal contributions, respectively. The Gaussian regulator is used. 





M [MeV] 


^J'I=l 


300 


350 


400 


450 


600 


val 


L 

NL 


2.910 
0.097 


2.519 
0.212 


2.339 
0.319 


2.245 
0.420 


2.174 
0.673 


sea 


L 

NL 


0.293 
0.122 


0.379 
0.198 


0.386 
0.238 


0.372 
0.262 


0.305 
0.289 


total 




3.422 


3.307 


3.282 


3.299 


3.442 



Table 5 

Contributions to the isovector magnetic moment, in units of the nuclear magneton. 
L and NL denote the local and nonlocal contributions, respectively. The Gaussian 
regulator is used. 

In Table ^ we list different contributions to the squared isovector magnetic radius, 
evaluated with the expressions given in App. |E.3| . This quantity depends on the 
prescription for the Noether current, but the difference between the straight- line 
path method and the weak-nonlocality approximation is even smaller than in 
the case of the baryon radius (see Table ^. The sea contribution is substantially 
larger than in the isoscalar electric case. The numbers are much larger than the 
experimental value of 0.77 fm^. 

To summarize the results for the observables, we note that the soliton, at the 
present mean-field treatment, is too large. This leaves room for the inclusion of 
other effects. One important effect comes from the center-of mass corrections 
which reduce both the mass and the mean square radii. The fields which describe 
the soliton break translational symmetry. The center of mass of the system is 
not at rest and it makes a spurious contribution both to the energy and to the 
mean square radii (more generally, to the form factors). This spurious contri- 



25 





M [MeV] 


{r^)Mj=i 


300 


350 


400 


450 


600 


val 


L 

NL 


2.498 
0.073 


1.288 
0.112 


1.043 
0.153 


1.001 
0.196 


1.285 
0.334 


sea 


L 

NL 


0.378 
0.137 


0.424 
0.187 


0.421 
0.220 


0.405 
0.243 


0.343 
0.276 


total 


straight line 
weak L 


3.087 
3.087 


2.011 
2.011 


1.838 
1.837 


1.846 
1.844 


2.238 
2.226 



Table 6 

Various contributions to the isovector magnetic mean square radius. L and NL denote 

the local and nonlocal contributions, respectively. The total result is displayed for the 

straight path and in the weak-nonlocality approximation. The Gaussian regulator is 

used. 

bution should be subtracted from the calculated values. The subtraction is a 
next-to-leading-order correction in N^. A rough estimate can be obtained from 
an oscillator model. If Nc particles of mass m move in the Is state of a harmonic 
oscillator of frequency hu, the center of mass of the system is also in a Is state and 
it contributes jhu: = (P^) /2Ncm to the energy. Thus, we can correct the soliton 
energies by subtracting (P^) /2Esoi from the calculated energy. Furthermore, in 
the oscillator model, the center of mass contributes a fraction -^ of the mean 
square radius, such that one corrects the mean square radius by multiplying the 
calculated value by a factor equal to (l — ;^j. These corrections to the soliton 
mass and the isoscalar electric radius may bring the calculated values close to the 



experimental ones |^8[. To some extent, the too large soliton may also reflect the 
lack of confinement of the model, or some other omitted dynamical factor, e.g. 



vector interactions ||3^ , ^ . Although the chiral field is sufficiently strong to bind 
the quarks, additional forces may reduce the soliton size. 



The problem with too low values of qa and of the isovector magnetic moment may 
be resolved by rotational 1/A'"c corrections, as found in soliton calculations which 
use proper-time or Pauli-Villars regularizations [^,0. The rotational corrections 
produce a factor lying between 1 and {Nc + 2)/Nc, which brings the calculated 
values closer to the experimental data. This problem requires a further study. 



6.4 Calculations with chiral fields constrained to the chiral circle 



In chiral quark models the soliton properties have so far only been calculated with 
the sigma and pion fields constrained to lie on the chiral circle |T2|,|39|. Although 
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Fig. 8. The energy of the soliton [in MeV] (bold lines), the three free-space quark mass 
(solid line) and the valence contribution to the soliton energy (thin lines) plotted as 
functions of M. Solid lines correspond to the unconstrained calculation (same as in 
Fig. 2), dashed lines to the calculation with the chiral-circle constraint, dashed-dotted 
to the calculation using proper time regularization, also with the chiral-circle constraint. 
In the latter case the energy of three free quarks is marked by 3M. The parameters of 
all models are fitted to obtain f-„ = 93 MeV. 



no clear physical ground for such a constraint can be seen in the derivations 
of effective quark models, it is nonetheless interesting to study its effect in our 
model. Fig. 8 shows the soliton energy of the constrained calculation as a func- 
tion of M. The constrained energy is 150 MeV higher than the energy of our 
solution at smaller values of M and stays some 70 MeV higher even in the region 
where the unconstrained chiral fields come closer to the chiral circle. From Fig. 3 
it is clear that in the unconstrained calculation it is energetically favorable to 
increase more the magnitude of the sigma field than the pion field. Contrary to 
the unconstrained calculation where the soliton energy at the threshold continues 
smoothly from the curve representing the energy of three free quarks, the energy 
of the constrained soliton starts abruptly at M = 298 MeV and in a small region 
of M just above the threshold stays above the 3Mg curve. Such an energetic insta- 
bility, absent in the unconstrained nonlocal model, is in fact common to several 



chiral models ||75|j76[| ; in particular, in the local model with the proper time regu- 
larization one finds a self-consistent solution for 338 MeV < M < 412 MeV even 



27 





M [MeV] 




300 


350 


400 


450 


600 


energy 

[MeV] 


Cval 

total 


327 
1249 


287 
1309 


279 
1336 


276 

1354 


275 
1387 


V /baryon 

[fm^] 


valence 
total 


1.190 
1.201 


0.711 
0.725 


0.667 
0.684 


0.677 
0.694 


0.798 
0.813 


9A 


valence 
total 


1.150 
1.205 


0.978 
1.079 


0.982 
1.084 


1.003 
1.100 


1.087 
1.164 


fJ'l=i 


valence 
total 


2.444 
3.223 


2.232 
3.065 


2.259 
3.088 


2.334 
3.146 


2.647 
3.388 


{r^)M,l=i 
[fm2] 


valence 
total 


1.692 

2.437 


0.873 
1.621 


0.844 
1.599 


0.930 
1.681 


1.465 
2.173 



Table 7 

Various observables (see Tables 2-6) obtained in the calculation with the chiral-circle 
constraint. "Valence" stands for the sum of the local and nonlocal contributions, the 
radii are calculated with the straight line prescription. 

though the sohton is energetically not stable in this region (see the dashed-dotted 
lines in Fig. 8). 

The observables calculated in the constrained calculations are displayed in Ta- 
ble 7. The differences to the unconstrained results presented in Tables 3-6 can be 
explained by two effects: due to a deeper effective potential in (|^) the valence 
orbits shrinks thus yielding a smaller valence contribution to the magnetic mo- 
ment and the two radii. On the other hand, due to the stronger chiral fields, the 
sea contribution to all quantities becomes more important and partially cancels 
the decreased valence contribution. 



7 Conclusions 



In this paper we have shown that a soliton, i.e. a bound state of A'"c = 3 quarks, is 
formed in a chiral quark model with nonlocal regulators. We have demonstrated 
its energetic stability and investigated its basic properties. Moreover, we have 
developed a scheme for quantizing the baryon number of the soliton, as well 
as for calculating observables. The construction of Noether currents has been 
accomplished through the use of straight line path-ordered P-exponents. This 
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construction is general and applicable to any model with nonlocal separable four- 
fermion interactions. 

We have shown that the nonlocal regularization, which is somewhat more and yet 
not prohibitively complicated, has several attractive features compared to the chi- 
ral quark models which use local regularizations, such as the Pauli-Villars or the 
Schwinger proper-time method. We have found that the pion field is considerably 
reduced compared to the local models which require the chiral field to remain on 
the chiral circle. The soliton is found to have properties which make it suitable 
for the description of the nucleon and for application of further corrections, such 
as projection [ ^ , ^J50|J5T|J5^ , or the inclusion of l/N^ corrections ||53| , |5^ |. 



The authors wish to thank Enrique Ruiz Arriola, Michael Birse, Klaus Goeke, 
Maxim Polyakov, and Nikos Stefanis for many useful discussions and comments. 



A The gauged nonlocal model 

Consider the gauge transformation of the quark field, 

g(x)-^e-'^"'^''(")g(a;), (A.l) 

where A" = 1/Nc for the baryon current, A*^ = r"/2 for the isospin current, 
and A" = 75r'^/2 for the axial current. We gauge the nonlocal model by using 
the path-ordered-exponent method described in Refs. pD| , pT| , |77| . The method is 



based on the Wilson line integrals (P-exponents) defined as 

where V denotes a path ordering operator needed for non-abelian gauge trans- 
formations, y^A^ is the gauge field (in general non-abelian), and s parameterizes 
an arbitrary path from x to y. The operator Wa{x, y) is a functional of the gauge 
fields y4^ with the following key property: when the gauge field undergoes the 
gauge transformation 

D^, -^ e-*^"*"D^e^^'''^% (A.3) 
where D^ = d^ + iA^A^^, the operator VF^(x, y) transforms as 

WAix, y) ^ Wa+8^{x, y) = e-'^''^'^^^^^WA{x, i/)e*""'^"(^). (A.4) 
Consider the quark-loop term of the Euclidean action (^), 

/g = -Tr In /5(-i9^ 7^ + m + r$r). (A. 5) 
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The gauged action is constructed by making the substitutions 

^M -^ ^M' (^ 1^1 y) -^ (^ I^aI y) = Wa (x, y) {x \r\y) , 

and 



p^ 



-i\"-<j,a (x) Q(^^i\"-<i>a (x) _ 



In exphcit form the gauged action reads 

/, {A) = -Tr In {-{Si-^^d^ + iS-^^X^Aa^, + I3m + VAP^rA) • 



(A.6) 



(A.7) 



(Ai 



Note that, in general, ta and (3 do not commute. In the gauge transformation the 
action ([A .81) transforms to 



/,(A)^/,(A + 90) 



— Trln 



-iPi.d, + /37m ^"^a^ + ^"'/'aM + [AVb, am:: 



+/3m + rA+a0e-^^'^^'"/3$e^^"*"r^+a^ 
Using the property rA+d(f> = e~*'*'"''^''r^e*'^'"''^" we find that 

'-^l>.^^ + 7mAM.^ + e*^"<^— -*""<^" 



(A.9) 



Ig{A + d^) = -Ti\ne-'^^'''^p 



'me 



+ 



(A.IO) 



Equation ( |A.1(J| ) shows that, provided e '^"Z? 



?Tie 



-iA''0a 



Pm (the case of baryon 
or isovector symmetries), or in the chiral limit m = 0, the action is invariant in 
the gauge transformation, 

/(A + 90) = /(A), (A.ll) 

such that the corresponding Noether current, evaluated from the expression 

SI (A) 



Jm ^^) = 
is conserved: 



6A^{x) 



A=0 



d,3l {^) = 0. 



(A.12) 



(A.13) 



B Explicit construction of Noether currents 



The action, expanded to first order in the A field, has the form 
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/(^) = -Tr- 



M„^ + r\f3<l>r + rf3<l>r\ 



E 



iu + ex (uj) 



(B.l^ 



where 



(x \rA\ y) = {x\r\y) il+ij ds^X^Aaf, (s) + ... = (a; \r\y) + (x \r\\ y) + ...(B.2) 
We can now write, with help of a representation of the 5 function, 



X 



y) = ij ds.X^A,^ {s)=ij -i^ j ds,e'^^-'>''X^Aa, (z) . (B.3) 



The Noether current is obtained from Eq. (|A.12|) , which, according to (|B.1| , |B.3| ) 
gives a general expression, valid for any path in the P-exponent: 



{uJ,X^\z)P'Ji,Xa{z\uJ,X^) 



(B.4) 



d4q 



(^, A^|x) {x\r\y) {y\XaP^r\u, X^) + {u, X^\rP^Xa\x) {x\r\y) {y\uj, X^) 

iuj + ex {uj) 

where j^^ and j^^ are the local and nonlocal contributions, respectively. We 
can now check explicitly the current conservation. With help of the equations of 
motion, we find 



dA (^) = E 



m{uj,X^\z){z\[(3,Xa]\uj,X^) 



\^ u; - lex [UJ) 



[B.5) 



E 



{uj,X^\z){z\(3Xar^r\uj,X^) - {LJ,XJ(3r^rXa\z){z\uj,X^) 



^ u-iex {UJ} 



P, -NL ( \_\^ {u,X^\z){z\(3Xar^r\u;,XJ) - {uj,X^\(3r^rXa\z){z\uj,XJ) 

{uj, X^ \r\z)Xaii^{z)i^z\r\ lo, XJ) - {uj, X^\r\z) l3^{z)Xa{z\r\uj , X^) 



E 



UJ — iex {oj) 



In the second equation we have used the identity 



d . r d^q 
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d^q ^^iiz-yyg _ ^^iz^xyg^ ^ ^(^^^ _ ^^ _ ^(^ _ ^^y ^g^^. 



(2vr) 



4 



Combining the local and nonlocal pieces ( p.5|JB.6D (which are not separately 



conserved), and using the equations of motion for the S and P fields, Eqs. 
we find that 

This immediately leads to the conservation of baryon and isospin current, and, 
in the chiral limit of ?7i = 0, to the conservation of the axial current. Note that 
the conservation laws are independent of the chosen path in the P-exponent. 

Quantities involving space integrals of Noether currents (charges, qa) also do not 
depend on the path. Indeed, the integration over z in Eq. (|B.4|) leads to 

X X 

and subsequently, through the use of the identity 

iiy^,-x^){x\r\y) = {x\rf,\y) , (B.IO) 

to 

f , . . ^ f du ^ (A^ I iPXqj^ + r^XaP<l>r + /Jr^Agr^) | A^) 

J J 2m Y^ uj — ie\ [uo) 

Note that any reference to the choice of the path has disappeared. Thus, charges. 



which are obtained from Eq. ( B.ll ) with /x = 0, or qa-, which has jjl = ?> {cf. Sect. 
[^ ) are independent of the path, hence are uniquely defined. One may easily 
generalize this result to any Green's function in the soliton background for the 
case where the external legs corresponding to Noether currents have vanishing 
four-momenta. The proof is straightforward through the use of the identity 



C Straight-line paths 



The expression ( p.4| ) is not, in general, suitable for calculating observables, since 
the path is not specified. A popular choice of the path is a straight line [ P(]| , pi| , 



Sf, = x^ + a{y^-Xf,) ds^ = da{y^-x^). (C.l) 
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Then in Eq. (|B.4|) we have 

/ (^ / •^''^''"""' = / (Ifj / ''" f^" - ""* <f'"-"-'''- (C.2) 

^ia(y,-.,)5(.-.-a(y-.)). (C.3) 



As a result, we find 







^.NL,straigh (^) 



± 

— ^2 '^4^ / '^4^ / '^'^ ^(^ — X — a{y — x)) X (C.4) 

(^,A^|x)(x|r^|y)(l/|Aa/j'l>r|u;, A,^) + {u, \JrP^\a\x){x\rf,\y){y\uj, X^) 

iuj + ca (cu) 

Since the fields $ are stationary, we can rewrite Eq. (p.4| ) in a simpler form 

j^t'''^'^''iz) = -I^Y.Id3xld,yJda5iz-x-a{y-x))x (C.5) 

A^ 

(A^|x)(f|r^(cu)|y)(y|A,/3$r(c^)|A^) + (A^|r(cu)/3'l'AJf)(x|r^(cu)|y)(y|A^) 

iuj + Ca (cj) 

In the second term above, we can interchange x with y and change the integration 
variable a ^ 1 — a. This leads to a manifestly Hermitian form. 



■NL.straigh /z* ' '^'^ 



- / — ^ / c/sx / d^y / da6{z — x — a{y — x)) x 

A^ 

(A^|x)(f|r^(^)|y)(y|A,/?'l 'rH|A^) ^ ^^ 



ZCU + Ca (Co") 

which is used below. 

D Weak-nonlocality approximation 



For the case where the nonlocality scale A is much larger than other scales in 
the problem {e.g. the inverse soliton size) we can commute the r^ operator with 
\x){x\ [0 in Eq. ( |C.6| ), thereby obtaining 
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■NL,weak / 
Jfjia \^ 



— ^ / dii^x / diy / da5{z — x — a {y — x)) 

+ h.c. 



{uj, \^\r^\x) {x\y) {y\\aP^r\uj, A^) 



tu + ex [uj) 

_ ^ {lu, X^\r^\z){z\\a/3^r\uj, \^) ^^ 

^ iu + exiuj) '^' 

Finally, commuting r^ and \z){z\ and using the stationarity of $ we find 

■NL,wcak (^\ = _ f^Y^ {\Jz){z\r^{u;)\ap^r{uj)\\^) 
•^^" ^ ' J 2-nXl iuj + ex{oo) ' ' ^ ' ' 

which is our current in the weak-nonlocality approximation. Note that some ar- 
bitrariness is involved here. We could have placed the \z){z\ operator anywhere 
between (A^| and |A^), and that would lead to different currents. However, all of 
these definitions become equal if the scale of the nonlocality is much larger than 
other scales in the problem, i.e. in the weak-nonlocality limit. 



E Evaluation of observables 

We calculate observables both with the straight-line prescription and in the weak- 
nonlocality limit. We work in the Minkowski space, by means of the replacement 

7o -^ i/3, joaix) -^ ipaix), (E.l) 

with Pa denoting the Minkowski charge density. 

E.l Baryon mean square radius 

The baryon charge mean squared radius equals to 

{r')^ = Jd,zz'pniz), (E.2) 

The contribution from the local part of the current is 
1 r dm f J\^\z)z^z\X^) 



B Nc J 2Tci T^ J UJ — iex [uj) 
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while the nonlocal term gives 



2\NL,straight 2 /", f duJ .^ [ ^ /", /", 

A^ Q 



^(i'-f-a(y-f))^^-'"^^">°(^^'^^<f^'^^"("^'^-\ (E.4) 

cu — ie\ [uo) 

We can now carry the a integration 

1 1 

d-iZ z^ / da5{z — x — a{y — x))= / da {x + a {y — x)Y = 



x^ + x- iy- x) + - iy- xf = -x^ - - {y - xf + -f (E.5) 

6 2 2 

to obtain 

„ \ NL, straight 111 , , 

r^) =-A + -B + -C, E.6 

/b 2 6 2' ^ ^ 

2 f duj ^ f ^ _^ {\Jro{uj)\x){x\P<^riu)\\^) 



iNr J 2-Ki^ J ^ 



^^ 2_ /"^v-y- {X^\rojj{uj)(3(!>r{uj)\X^) 



iN^ J 27ri '^p{ u - iex [uj) 

2 f du ^^ f 2(^c.|^)(^ko(^)/3'^^(^)|Ac,) 



iNc J 27ii Y^ J ui — iex {^) 

In the derivation of B we have used the identity 

3 

{x - yf {x\ro{u;)\y) = -Y,{x\rojj{uj)\y). (E.7) 

i=i 

The above formulas can be further decomposed into the valence and sea contri- 
butions: 



(^^)^' =z^^ijd3zz^\{z\val)\'^, 

Nc J 2tt r^ u^ + ei (u) J 



r 

B 



A™ =4zvaievai / d^xx (val|r'(ievai)|f)(f|/3<^r(ievai)|val), 

^^^^ = 1^ /I^E 2f2i fdsxx'{X^\r'{uj)\x){x\PM^)\X^), 
Nc J 271 Y^ u^ + ex [uj) J 

E^^' = 16z.,ie.,i(val|(-V')r'"(2e.ai)/9$r(^e.ai)|val), 
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B 



16 r dui 



E 



UJ 



NcJ 2n^u}^ + ei{uj] 



(AJ(-V^)r'"H/3$r(.;)|A^), 



2 I „2 

C™' = 4zvaievai / (i3xx^(val|f)(f|r'(ievai)/?$r(ievai)|val), 

2 






ci3a;x2(A^|x)(x|r'H/3$r(a;)|A^).(E.8) 



In the derivation we have used the notation r'(fc^) = d/dk"^ r{k'^), etc., and the 
fact that r{uj) = r{—uj). 

In the weak-nonlocality approximation we obtain 



r^) 



NL,weak 



B 



C. 



(E.9) 



E.2 QA 



We evaluate qa from the expression 
2 



Z f 
9a = --J dszAssiz), 



(E.IO) 



where — ^ is a factor coming from the cranking projection method |^ and ^33 is 



the hedgehog matrix element of the third-space, third-isospin component of the 
axial charge. According to Eq. ( |B.11|) we find 



9A = - I 7^2^ — 7-^-^ , (E.ll) 



2m 



Xu. 



LU — tex [LUj 



which gives 



L,val 



Nr 



Qa = — ^Zvai ( val 1 0-3X3 1 val) , 
L,sea^l f duj^ ex (c^) (A^ [0-3X31 A^) 



NL,val 

9a 



NL,sea 



AN, 



Zvai(val (-iV3)r'75r3/5$r val) 



9a 



2tx 



E 



duj ^ ca [uj] {X^ {-iVsy-f^T^p^r A^) 



u;^ + el (a;) 



(E.12) 
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E.3 Isovector magnetic moment and mean square radius 



The isovector magnetic moment is obtained from the expression 

/i/=i = -le^'"" fdsz z,V3m{^), (E.13) 



where the — | is the cranking projection factor [57,7^] and V^m is the hedgehog 
matrix element of the m-space, third-isospin component of the isovector current. 
For the local contribution we find 

N r 
th=i^ = — T^z^ai / dsz {va\\z)(3{z x 7)3r3(ilval), (E.14) 

L,sea I f duJ^^ex{uj)Jd3z{XJz)l3{zX^)3Ts{z\X^) 

6 J ZTT Y^ uj^ + el [u) 

For the nonlocal contribution with the straight-line prescription we have 

3 m ^ 

^NL,straight ^ __^ I ^^^ ^^ I ^^^ I ^^^ / ^q, ^(^^ _ ^ _ ^ (^ _ ^)) ^^^ 



J 2m Y^ uj — ie\ [u) 

Then 

_3jrra 



^NL.straight ^ __^ I ^^^ I ^^^ I ^^(^^ j^ a {y - x))j Rm 



= -7- j dsx J d^yixj + - [jjj - Xj)) Rm- (E.17) 



Through the use of the identities {yj — Xj) {x\rm{uj)\y) = —i{x\rjnj{uj)\y) and 



^^■''^fmji^) = we obtain 

NL,straight ^^^"^ f, f duj ^^ {XJx)Xj{x\rm{uj)T3f3^r{uj)\X^) 

/ii=i = -^ d^x Tr-l^ ■ — r~\ + ^•^• 

(E.18) 

Since e^^"^Xjrm{uj) = e^^"^Xj{—2i'Vm)r'{u) = 2L3r'{uj), with L being the orbital 
angular momentum operator, we can write 

NL,straight _ 1 f duj ^^ {X^\L3rmiuj)T3f3^r{uj)\X^) 



m. straight ^t :^y \-'a.|^.>mV-r3P^ V^;K'W ^ ^ -^ 

^^-^ 3 J 2m Y u}-iex{uj) ^ ' 
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Now we notice that the same expression follows when we use the weak-nonlocality 
approximation, therefore 

NL, straight NL,wcak /t-\ on\ 

l^i^i = f^iJi (E.20) 

Finally, 

2N 
^NL,vai ^ -^Zvai (val|r'(?evai)-^3T3/3$|r | val) , (E.21) 

Since [r', L3] = 0, we are free to interchange the order of r' and L3 in the above 
formulas. 

The isovector magnetic mean square radius is defined as 

(r^)^^^ j^^ = -^/'"^ J d,z z,z^V^^[z), (E.23) 

For the local part we find 

(r^)"^'™' = -TT^Zvai fdszz^Yii\\z)P{zx 7)3r3(f|val), (E.24) 

\ /m,I=l 6/i/=i J 

/ 2\L,sea _ 1 f duj y. cx juj) J dsz z"^ {\^\z} p{z X ^)3r3{z\\^) cj^ 25) 
V /m,i=i ~ 6/i/=i ^ 27r ^ ct;2 + ei (cu) ' I • J 

For the nonlocal part we have 

/ d^z z^Zj / (^30; / d^y / (iQ;5(z — x — a {y — x)) Rm 

^"^ 

da{x + a {y — x))^(x + a {y — x))j Rm (E.26) 



„\NL,straight e^"^ 

r^\ = — 

I m,I=l 6Uf_ 

^ - Q 



6/i/=i p 



rfsa; / d^y 



\ 2 ila _l/'^_^^2 
r^x Xj -\- y Xj [y X) Xj 



R. 



m- 



Finally 

„\ NL.straight 111 , , 

r^)^^ j^^ = -Am + -Bm + -C^, (E.27) 

^™' = -^^z,ai / rf32;z2(val|r'(ie,ai)i^3|^r3/^<l>(?)(f|r|val), 

.sea^ 2 /•da;^ eAM/c?32:2:^(A^|rXa^)i^3|^T-3/3$(f)(f|r|A^) 
"^ 3/i,=ii 27rf- a;2 + ei(cu) 
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QAJ 

BZ' = -7r^Zvai(val|(-V')r'"(2e.ai)i:3r3/3$r|val), 
o/i/=i 

Rsea^^ (du^ e^ juj) {XJ{-Vy"{uj)L,TsP<!>r\X^) 

™ Q,,. / O^ 2^ , ,2 I „2 



3/i/=i J 27r Y^ cu^ + e| (cj) 

In the weak-nonlocality limit 

„\ NLjWeak 
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